Quark matter under strong magnetic fields in the su(3) Nambu— Jona-Lasinio Model 
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In the present work we use the mean field approximation to investigate quark matter described 
by the su(3) Nambu-Jona-Lasinio model subject to a strong magnetic field. We consider two cases: 
pure quark matter and quark matter in /3-equilibrium possibly present in magnetars. The results 
are compared with the ones obtained with the su(2) version of the model. The energy per baryon 
of magnetized quark matter becomes more bound than nuclear matter made of iron nuclei, for B 
around 2 x 10 19 G. When the su(3) NJL model is applied to stellar matter, the maximum mass 
configurations are always above 1.45 Mq and may be as high as 1.86 Mq for a central magnetic 
field of 5 x 10 18 G. These numbers are within the masses of observed neutron stars. 
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I. INTRODUCTION 

In non-central heavy ion collisions such as the ones per- 
formed at RHIC and LHC-CERN, physicists have been 
looking for a possible signature of the presence of CP-odd 
domains in the presumably formed quark-gluon plasma 
phase The study of deconfined quark matter sub- 
ject to strong external magnetic fields is then mandatory 
if one intends to understand the physics taking place in 
such colliders. 

Neutron stars with very strong magnetic fields of the 
order of 10 14 — 10 15 G are known as magnetars and they 
are believed to be the sources of the intense gamma and 
X rays detected in 1979 @,||. The hypothesis that some 
neutron stars are constituted by unbound quark matter 
cannot be completely ruled out 0] since the Bodmer- 
Witten conjecture @ cannot be tested on earthly experi- 
ments. This conjecture implies that the true ground state 
of all matter is (unbound) quark matter because theoret- 
ical predictions show that its energy per baryon at zero 
pressure is lower than 56 Fe binding energy. 

In the present work our aim is to investigate quark 
matter described by the su(3) version of the Nambu- 
Jona-Lasinio ((J model exposed to strong magnetic fields. 
In the case of pure quark matter, as predicted by the 
QCD phase transition possibly taking place in heavy ion 
collisions, the magnetic field is certainly external. In the 
case of neutron stars, the magnetic field can be generated 
by the alignment of charged particles that are spinning 
very rapidly. We next use an external field to mimic the 
real situation, which we do not know how to determine. 
Albeit in an approximate way, the effect of the magnetic 
field on the macroscopic quantities as radius and masses 
can be obtained. 

Recently the su(2) version of the NJL model was used 
to treat both situations described above We have 
shown that, for pure quark matter, the energy per baryon 
for magnetized quark matter has a minimum which is 
lower than the one determined for magnetic free quark 



matter. We have also obtained that a magnetic field of 
the order of 2 x 10 18 G barely affects the effective mass 
as compared with the results for matter not subject to 
the magnetic field. For B = 5 x 10 19 G matter is to- 
tally polarized for chemical potentials below 490 MeV. 
For small values of the magnetic fields the number of 
filled Landau levels (LL) is large and the quantisation 
effects are washed out, while for large magnetic fields the 
chiral symmetry restoration occurs for smaller values of 
the chemical potentials. When /^-equilibrium is enforced, 
the numerical results show that, for the the su(2) case, 
only very high magnetic fields (B > 10 18 G) affect the 
equation of state (EOS) in a noticeable way. 

The inclusion of the s-quarks, necessary in the su(3) 
NJL model, poses some new numerical difficulties and 
some questions that need to be addressed. Those prob- 
lems are tackled through out the paper. One of the ques- 
tions was raised in [J Q and refers to the stability of 
quark matter described by the NJL model. The authors 
show that it is not absolute stable. As already mentioned, 
in we have seen that the inclusion of the magnetic field 
increases stability in the su(2) version and the same be- 
havior is expected in the su(3) NJL, which is shown next. 

The paper is organized is such a way that all calcula- 
tions already shown explicitly in 0] are not repeated but 
all important differences arc outlined. In sections II and 
III the formalism (mean field theory) and the equations 
of state are shown and in section IV the final results are 
displayed and the conclusions are drawn. 



II. GENERAL FORMALISM 

In order to consider (three flavor) quark stars in f3 
equilibrium with strong magnetic fields one may define 
the following lagrangian density 



(1) 



where the quark sector is described by the su(3) version 
of the Nambu-Jona-Lasinio model 

£f = $f hn - q f A») - rh c ] yj f + C sym + C det , (2) 



where C sym and Cdet are given by: 

8 

Csym = G^[(^ / A a ^ / ) 2 + % i 75A Q V/) 2 ] , (3) 

C de t = -K {detj [^/(l + 76)^/] +det y [^/(l - 75)^/] } 

where if>f = (u, d, s) T represents a quark field with three 
flavors, to c = diagy (m u , m,d, m s ) is the corresponding 
(current) mass matrix while qf represents the quark elec- 
tric charge and A denotes the Gell-Mann matrices. Here, 
we consider m u = m,d 7^ Tn s . The Cdet term is the t'Hooft 
interaction which represents a determinant in flavor s pac e 
which, for three flavor, gives a six-point interaction [ll| 

det/$/C%) := J2 e t3k (uO^)(dOiP 3 ){s0^ k ) , (5) 

where e^-fe is the usual three-dimensional Levi-Civita 
symbol. The lagrangian also contains the C sym term 
which is symmetric under global U{Nf)i, x U (Nf)^ trans- 
formations and corresponds to a 4-point interaction in 
flavor space. In the appendix we discuss the steps to 
obtain Cf in the mean-field approximation (MFA). 
The leptonic sector is given by 

Ci = to [ 7ju (id" - qiA") - mi] to , (6) 

where I = e, /z. One recognizes this sector as being repre- 
sented by the usual QED type of lagrangian density. As 
usual, and i*) t „ = d^A v — d v A^ are used to account 
for the external magnetic field. Then, since we are in- 
terested in a static and constant magnetic field in the z 
direction, A^ = S^xiB. 

III. THE EOS 

We need to evaluate the thermodynamical potential 
for the three flavor quark sector, Qf, which as usual can 
be written as Qf = —Pf = £f — TS — HfPf where 
Pf represents the pressure, £f the energy density, T the 
temperature, <S the entropy density, and \i f the chemical 
potential. 

For the present study, just the zero temperature case 
is important and, as a consequence, the term with the 
entropy vanishes. The total pressure for three flavor in (3 
equilibrium is given by 

P(Hf,H h B) = Pf \ Mf + If\ mi + !L , (7) 

where our notation means that Pf is evaluated in terms 
of the quark effective mass, Mf, which is determined in 
a (nonperturbative) self consistent way while P; is eval- 
uated at the leptonic bare mass, m;. The term P 2 /2 
arises due to the electromagnetic term P Mt/ P A " y /4 in the 
original lagrangian density. The subscript N indicates 
normalized pressures. Here, our normalization choice is 
such that Pf = Q at /// = (/ = it, s, d) and P t N = at 
Hi = (l = e, n) implying that P(0, 0, B) = B 2 /2. 
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A. Quark Contribution to the EOS 

In the mean field approximation the pressure can be 
written as 

Pf = e u + e d + e s - 2G(^ + $ + ft) + ak^^ , (8) 

where an irrelevant term has been discarded. The pres- 
sure due to the three quarks is diagrammatically repre- 
sented in figure la. 
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(b) 

FIG. 1: a) Feynman diagrams contributing to the quark pres- 
sure in the MFA. The lines represent the three dressed quark 
propagators for the different flavors: u (continuous line), d 
(dashed line) and s (dotted line) . The black dot represents G 
and the black hexagon represents K. b) Diagrammatic rep- 
resentation of the effective mass for flavor u. The diagrams 
contributing to the other two flavor display the same topology. 

For a given flavor, the Of term is given by 
and the condensates, <j>f are given by 

»/-<*/*>— '/g& V: O.0 (10) 

where all the traces are to be taken over color (N c = 3) 
and Dirac space, but not flavor. In order to obtain results 
valid at finite T and [L in the presence of an external 
magnetic field B one can use the following replacements 

po — > i(u v - ifJ-f) , 
p 2 -> p 2 z + (2n + 1 - s) , with s = ±1 , n = 0, 1, . . . 




In the above relations, uj v = (2v + 1)ttT, with v = 
0, ±1, ±2, . . . representing the Matsubara frequencies for 
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fermions while n represents the Landau levels (LL) and 
s represents the spin states which, at B ^ 0, must be 
treated separately. The case T = in which we are inter- 
ested can be easily obtained after the above substituions 
(seeRef. 0). 

The effective quark masses can be obtained self consis- 
tently from (see figure lb) 



Mi 



(11) 



^ = (0p c + ^7 a9 + ^ ed ) A// 



where 



h vac _ MfN c 



Ae A - Mj In 



A + e A 



A I 



f 



(17) 



, (18) 



with (i,j,k) being any permutation of (u,d,s). So, to 
determine the EOS for the su(3) NJL at finite density and 
in the presence of a magnetic field we need to know the 
condensates, <pf, as well as the contribution from the gas 
of quasiparticles, Of. Both quantities, which are related 
by (j>f ~ d9f/dMf, have been evaluated with great detail 
in Ref. [7j- Here, we just quote the results 



_ ^ py ac _|_ p ma 9 _|_ pmeil 



where the vacuum contribution reads 



(12) 



pvac = _^) M j ln 
/ 87T 2 ' 



(A + e A ) 



- £A A(A 2 + e |) 



(13) 

where we have defined e A = \j A 2 + M 2 with A repre- 
senting a non covariant ultra violet cut off. The evalua- 
tions performed in Ref. 7;] also give the following finite 
magnetic contribution 



, mag = N c {\q f \Bf 
f 2tt 2 



C\-l,x } ) - \{x) - x f )\nx f + 



(14)" 

where x f = Mf/(2\q f \B) while ('(~l,x f ) = 
d((z,Xf)/dz\ z= -i where Xf) is the Riemann-Hurwitz 
zeta function (lp| . Finally, after integration, the medium 
contribution can be written as 



Tjmed 



E 

k=0 



Qk 



\df\BN c r 



4tt 2 



s f (k,B) 2 hi 



H f Jfij - s f (k,B) 



A 4 / ~ s f{k, B) 2 



s f (k,B) 



as) 



where Sf(k,B) = ^ M 2 + 2\qf\Bk, a = 1, a k> o = 2. 
The upper Landau level (or the nearest integer) is defined 
by 



Pi 



M 



r2 



Pj,f 

2\q f \B 2\q f \B- 



(16) 



Finally, the condensates 4>t entering the quark pres- 
sure at finite density and in the presence of an external 
magnetic field can also be written as 



lUIUlj 



M f \q f \BN c 
2tt 2 



[\nT(x f ) 



±\n(2ir) + x f - ±(2x f - l)\n(x f ) 



,(19) 



and 



ijned 



E 

fe=0 



ln 



Oik 



M f \q f \BN c 
2tt 2 



H f + ^Jli 2 f -Sf{k,Bf 

8f(k,B) 



(20) 



From the pressure one can obtain the density, pf, corre- 
sponding to each different flavor, which is given by 



Pf 



\q f \BN C u 



k=0 



(21) 



where kpj = yM/ — s f{k,B) 2 , since dP/d(j>f = 0. 
The quark contribution to the energy density is 



£ f (fi f ,B) = -P? + J2»fPf - ( 22 ) 



where Pf = Pf{^f)\ Mf {^) ~ P/(0)| M/ (o)- 

Throughout this paper we consider the following set of 
parameters [HI: A = 631.4 MeV , m u = nid = 5.5MeV, 
m s = 105.66MeV, GA 2 = 1.835 and KA 5 = 9.29. 



B. Lepton Contribution to the EOS 

The leptonic contribution, P; has also been evaluated 
in detail in Ref. Q where the normalization requirement 
Pj W = at fii = has been adopted. The result shows 
that, at the one loop level, only the following (finite) 
medium contribution has to be considered 



^ = E E 



Oik- 



l=e fe=0 



4tt 2 



PiVPi - si(k,B) 2 



si(k,BY In 



2 , | Mi + VM? - si{k,B) 2 



si(k,B) 



(23) 
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Then, the leptonic density is also easily evaluated yield- 
ing 



pi= 2_ ot k —^-k F .i(k,Si) , 



fc=0 



(24) 



where kpj{k, s{) = y/x? — si(k, B) 2 . Finally, the leptonic 
energy density reads 



£ l (p. l ,B)=-P? + Y l HiPi 
i 



(25) 



The lepton masses are m e — 0.511 Me V and m M 
105.66 MeV. 



IV. RESULTS AND CONCLUSIONS 

In the sequel we consider two different situations of 
quark matter under a strong magnetic field: a) pure 
quark matter with the same chemical potential for all 
quark flavors; b) /3-equilibrium quark stellar matter. 

We first discuss the properties of pure quark matter 
with equal chemical potentials for all flavors, namely 
the behavior of the dynamical quark masses, the chiral 
symmetry restoration with density and the energy per 
baryon. In Fig. [2] we display the masses of quarks u 
and d as function of the chemical potential for different 
values of the magnetic field and the two versions of the 
NJL model. For the magnetic field intensities used, one 
can clearly identify the filling of different Landau levels 
causing the usual kinks in the curves. For the three in- 
tensities considered the chiral symmetry is approximately 
restored for fi = 400 MeV. 

It is interesting to see that although the general behav- 
ior is the same, the effect of the LL is more pronounced 
in the su(2) version. 

In Fig. [3] the mass of the s quark is shown as a func- 
tion of the chemical potential for different values of the 
magnetic field. One can see how drastically it falls around 
H = 450 MeV. For magnetic free quark matter, this is the 
same behavior shown in Fig. 3 of (T^ |. One can observe 
that the curve is no longer smooth when B is turned 
on, but the values of the strange quark mass do not 
vary much. According to [13], the fact that the strange 
quark mass remains relatively high as compared with the 
masses of the other two quarks is the main reason why 
deconfined quark matter may not be likely to appear in 
the core of hybrid neutron stars. For a magnetic field 
larger than 10 19 the restoration of chiral symmetry for 
the s-quark occurs in steps and starts at a smaller chem- 
ical potential than the B=0 case. 

The phenomenon of magnetic catalysis, which en- 
hances chiral symmetry breaking, has been well discussed 
within the su(2) version of the NJL model [3]. Here, for 
reference, we show the vacuum effective mass of the three 
quarks as a function of the magnetic field in Fig. 3J For 
B > I0 19 G the vacuum masses increase dramatically 




320 340 360 380 400 420 440 460 480 

H(MeV) 

FIG. 2: Mass of the quarks (a) u and (b) d as a function of 
the chemical potential for B — 0, 
the su(2) NJL and su(3) NJL. 



10 iM and 2 x 10 1J G within 



with the magnetic field as expected. A similar increase 
of the vacuum mass was also obtained for the su(2) ver- 
sion of NJL in 0, EH and the effect is related to the fact 
that the B field facilitates the binding by antialigning 
the hclicities of the quark and the antiquark, which are 
then bound by the NJL interaction. As shown in Fig. 
31 an interesting result of the su(3) version is that, due 
to its larger electric charge, the u quark has an effective 
mass that becomes larger than that of the s quark for 
B > 1.5 x 10 20 G. 
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FIG. 3: Mass of the s quarks as a function of the chemical 
potential for B = 0, 10 19 and 2 x 10 19 G within au(3) NJL. 

In Fig [5] the baryonic density is shown as a function of 
the quark chemical potential for two values of the mag- 
netic field and for both versions of the NJL model, su(2) 
and su(3). As already noticed in 0], once again, for small 
values of the magnetic fields the number of filled LL is 
quite large and the effects of the quantization are less 
visible. Due to the Landau quantization, the increase of 
the strength of the magnetic field gives rise to a decrease 
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FIG. 4: Vacuum mass of the quarks as a function of the mag- FIQ g . Energy pgr nudeon ag & function of density for B = 
netic field B. 0) 1Q i9 and 2 x 1Q i9 G within NJL su ( 2 ) and NJL su(3). 



of the number of the filled LL and the amplitude of the 
oscillations is more clear in the graphics. For each value 
of the magnetic field, the kink appearing at the smallest 
chemical potential corresponds to the case when only the 
first LL has been occupied. 
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FIG. 5: Baryonic density as a function of the quark chemical 



potential for B = 0. 
and su(3) NJL. 



10 1J and 2 x 10 19 G within both su(2) 



In Fig. [5] one can see that the inclusion of the mag- 
netic field makes matter more and more bound in both 
versions of the model. For the present set of parameters, 
the energy per baryon E/A of magnetized quark matter 
becomes more bound than nuclear matter made of iron 
nuclei, f \™ Fe ~ 930 MeV for B around 2 x 10 19 G. 

We next consider stellar matter made out of quarks, 
electrons and muons in /^-equilibrium, as possibly occur- 
ring in the interior of magnetars. It is worth mentioning 
that, in this case, the three different quarks bear different 
chemical potentials, determined by the chemical equilib- 
rium conditions 

Pd = Ps = Pu + Me, Mm = Me- 

We start by plotting the quark effective masses for dif- 
ferent values of the magnetic field in Fig. \7\ It is seen 
that the results for non-magnetized matter (B = 0) al- 
most coincide with the ones obtained for B — 10 18 G. A 
decrease of the s quark mass starts only at ~ 0.8 mi -3 . 
This behavior had already been discussed in [l3| . If the 



magnetic field is strong enough the mass of quark s oc- 
curs in finite jumps which may give rise to an increase of 
the strangeness fraction as shown in Fig. [5J 
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FIG. 7: The quark effective mass for /3-equilibrium quark 
matter with a constant magnetic field within NJL su(3). 

The quark fractions Yi = pi/p, i = u,d,s are shown 
in Fig. [5J Again the results for B = are similar to 
the ones for B = 10 18 G. For strong enough fields the 
quark u fractions increase with a reduction of the quark 
d fraction. The quark s fraction has a sudden increase 
for p ~ 0.7 fm~ 3 but above p ~ 0.9 fm~ 3 remains below 
the B = fraction. 
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FIG. 8: Fraction of quarks in /3-equilibrium quark matter for 
a constant magnetic field within NJL su(3). 
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In Fig. [9] the EOS for different values of the magnetic 
field is shown. For magnetic fields as large as B — 10 18 G 
the differences are very small as compared with non- 
magnetized matter. For larger fields there is an overall 
net softening of the EOS. 
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FIG. 9: Equation of state of /3-equilibrium quark matter for 
constant magnetic fields within the NJL su(3). 

It is well known that at the surface the magnetic field 
should not be larger than ~ 10 15 G. We have introduced 
a density dependent magnetic field as in @, HH : 



Bi = B surf + B Q 



where B SXLI ^ — 10 15 G is the magnetic field at the sur- 
face, Bi is the magnetic field at the interior of the star 
for large densities and the parameters a — 5 x 10~ 5 and 
7 = 3 were chosen in such a way that the field increases 
fast with density to its central value but still describes 
correctly the surface of the star where the pressure is 
zero. We show the equations of state for quark matter 
in /3-equilibrium and a density dependent magnetic field 
within both versions of the NJL model in Fig. [TDJ As 
implicit in eq. (|26|) . the field at the surface is 10 15 G. 
The magnetic field makes the EOS harder with conse- 
quences in the gravitational and baryonic masses of com- 
pact stars, whose properties are obtained from the in- 
tegration of the Tolman-Oppenheimer-Volkoff equations, 
which use as input the EOS obtained with the density de- 
pendent magnetic field. The results are displayed both in 
Fig. [TTIand in TablclH from where it is seen that both the 
gravitational and the baryonic masses increase with the 
increase of the magnetic field for an intensity larger than 
- 5 x 10 18 G for the su(3) version and 10 18 G for the su(2) 
NJL. However, the increase of the gravitational mass is 
larger than the increase of the baryonic mass because 
the contribution of the magnetic field becomes more and 
more important as the field increases. This explains the 
decrease of the central energy/baryonic density for the 
stronger fields considered. 

Another important effect of the field on the properties 
of the stars is the increase of the radius of the star with 
the largest radius, which may be as high as 9.5 Km for 
the su(3) NJL. In general, the maximum mass star con- 
figurations for the su(2) version of the NJL model are 



1 — exp 



Po 



(26) 



smaller with smaller radius, ~ 7 Km, in average 2 Km 
smaller than the corresponding stars in the sit(3) version 
of the NJL model. 

Within the su(3) NJL the maximum mass configura- 
tions are always above 1.45 M@ and may be as high as 
1.86 Mq for a central magnetic field of 5 x 10 18 G. These 
numbers are within the masses of observed neutron stars. 
On the other hand the su(2) version of the NJL model 
forsees too small star masses except for very large mag- 
netic fields. 




FIG. 10: Equation of state of /3-equilibrium quark matter for 
a density dependent magnetic field within NJL su(2) and NJL 
sw(3). The EOS for B = is also shown. 
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FIG. 11: Mass-radius curves for the families of stars within 
NJL su(2) and NJL su(3) corresponding to the EOS shown 
in Fig. [9] For B < 10 18 the curves coincide with the B — 
results. 

The effects of the anomalous magnetic moments has 
been shown to be relevant Il7l . Il8l | and we intend to 
take them into account in the next calculations. 

The color superconductivity (CS) which allows 
the quarks near the Fermi surface to form Cooper pairs 
that condense and break the color gauge symmetry (20| 
is known to be present in the QCD phase diagram at 
sufficiently high densities. The effect of strong magnetic 
fields on the CS properties of quark matter, which can 
be drastic for sufficiently high fields, has already been 
studied by several authors [2l| . It would be important to 
investigate how this SC phase could affect the properties 
of quark stars under strong magnetic fields. However, 
it could be that CS is only affected by magnetic fields 
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TABLE I: Maximum mass configurations for NJL sw(3) and 
su(2), and several magnetic field intensities: the gravitational 
mass (M), baryonic mass (Mb), radius (R), central energy 
density (e c ), baryonic density (p c ) and magnetic field (B c ) 
are given 
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1.19 


10 15 


10 18 


1.46 


1.53 


8.93 


7.49 


1.19 


1.6xl0 17 


5xl0 18 


1.47 


1.54 


8.88 


7.94 


1.24 


8.8xl0 17 
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1.38 
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1.72 


4.0xl0 18 
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1.49 


1.41 
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14.47 


1.49 


5.7xl0 18 



where tr is the trace operator in flavor space. So, adding 
and subtracting these expressions, we can rewrite the 
NJL symmetric four-point interaction term as: 

8 

C sym = G^[(lpf\ a 1pf) 2 + (lpfij 5 X a ^f) 2 ] 
8 

= tr(\ a $) tr(\ a &) = 2G tr(<S> $ f ) . (A2) 

a=0 

The summation involved in the latter equality can be per- 
formed noting that an arbitrary matrix A in the Nf=3 
flavor space, can be expanded in terms of Gell-Mann ma- 
trices as follows: 

8 1 
A = X/ CaXa ' with Ca = 2 tr ^ XaA) ' ^ A3 ^ 

a=0 

The expansion coefficients c are obtained using the the 
Gell-Mann matrices property: tr(\ a \ b ) — 26 a b- So, we 
can write: 



stronger than the ones considered in the present paper, 
which, however, predicts already a very high maximum 
mass, M ~ 1.9 M@. The largest magnetic field we got 
in the center of a quark star is 5 x 10 18 G, while in 21 1 
it is shown that a noticeable effect requires fields above 
~ 10 19 G. 
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APPENDIX A: THE SU(3) NJL MODEL IN THE 
MFA 



In this appendix the main steps in order to obtain the 
Nambu- Jona-Lasinio lagrangian given eq. ^ in the mean 
field approximation are explicitly shown. Firstly, we 
consider the C sym term given in eq.([3]). For later conve- 
nience, we define the matrix elements of $ and its adjoint 
$t as [2|: 

$y = ipjO- - 75 M . ®\j = + 75)^ , 

where i, j are flavor labels. From these definitions, one 
can easily show that: 



t/jf(l -7s)A a V'/ = tr(X a ^) , 
^(1 + 75)A„V/ = tr(X a ^) , 



(Al) 



8 8 1 8 

tr(AA^) =tr(J2caKj2 c t X h = ^E'^^'^^ 4 ^ 



a=0 



b=0 



a=0 



(A4) 

where in the latter term we have used that the Gell-Mann 
matrices are hermitian, i.e., A a =AJ. We then evaluate 
C sym m the mean field approximation linearizing the in- 
teraction terms. We follow refs.Jll], [13] approximating 
the product of two operators 0\ and O2 by: 

Oi6 2 «Oi(0 2 ) + (Oi)0 2 - (6i>(d 2 ) ■ (A5) 

Therefore, calculating explicitly the trace involved in 
eq. (|A2[) and taking into account the prescription above, 
C sym , can be written in the MFA as: 



C 



syrn 



AG 



<p u v)u + <pd cftd + 4> s s^i 



1 



where we have used: 

(ipitpj) = Sijfa and (ipa5^j) = 



(A6) 



(A7) 



The only 3 non-vanishing terms are the condensates 
which were defined in eq. (|10|) . Finally, we consider the 
t'Hooft term, eq.Jl]), which is a six-point interaction in 
the su(3) flavor space. Notice that term involves the 
product of three operators which we linearize analogously 
to ea.([A"5)l: 



(Oi)0 2 (0 3 ) + (Oi)(0 2 )0 3 



O1O2O3 wOi(0 2 )(0 3 ) 
- 2(0 1 )(6 2 )(0 3 ) . 



So, in the MFA the determinants which appear in the 
t'Hooft term can be written as: 



8 



Y,eijk[fr04i)(dOil>j)(aOil> k ) + (uO^)(d0^ 3 )(s0^ k ) 

+ (uO^)(dO^)(s0^ k ) - 2(uO^)(dO^)(s0^ k ) . 

Now, inserting the operator O — 1 ± 75 and using the 
properties given in eq. (|A7|) . we obtain the t'Hooft term 
in the MFA : 

C-det = -2K{(j) d (f) s uu + <j) u <j) s dd + 4> u (t)dss - 2<p u (j>d4 l s) ■ 

(A8) 



From eq.© and eqs. (|A6IA8[) the su(3) NJL la grangian 
in the MFA is given by: 

cf FA = $ f ( 7 „ w - if An - m) ti>f 

-2G{4>1 + <& + </%) + iKfaMs , 



where M is a diagonal matrix with elements defined in 
eq,(fTT|). 
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